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INTRODUCTION 

The  active  control  of  fully  developed  turbulent  flows  are 
of  particular  interest  for  many  industrial  applications.  In 
fact,  as  regards  aerodynamic  applications,  the  control  and 
modification  of  turbulent  flows  could  be  used  either  in  ex¬ 
ternal  flows,  to  reduce  the  total  drag  of  aerospace  vehicules 
or  in  internal  flows,  to  increase  the  mixing  rate  in  a  scram- 
jet  for  instance.  However,  as  these  industrial  applications 
are  complex,  we  limit  the  present  study  to  the  analysis  of 
a  basic  free  turbulent  flow:  the  incompressible,  turbulent 
mixing  layer  developing  spatially  downstream  a  flat  plate. 

In  such  a  flow  field,  large-scale  coherent  sructures  exist 
which  contain  most  of  the  turbulent  kinetic  energy.  These 
coherent  structures  are  mainly  responsible  for  vibrations, 
noise  generation,  etc...,  Therefore,  in  term  of  control,  it 
seems  important  to  describe  correctly  the  characteristics 
of  these  structures  and  to  predict  precisely  their  time  evo¬ 
lution  using  models  as  simple  as  possible.  One  of  the 
methods  proposed  to  mimic  the  dynamical  behaviour  of 
the  flow  is  to  develop  a  low-order  dynamical  system.  This 
study  presents  the  ability  of  the  Proper  Othogonal  Decom¬ 
position  (POD)  (Lumley,  1967)  to  obtain  this  low-order 
dynamical  system.  In  the  POD  approach,  the  coherent 
structures  are  defined  in  terms  of  optimal  signature  of  the 
turbulent  kinetic  energy.  The  coherent  structures  are  then 
associated  to  the  eigenfunctions  of  the  POD.  This  method 
has  already  been  used,  with  success,  to  study  either  the 
near-wall  evolution  of  the  flow  within  a  turbulent  bound¬ 
ary  layer  (Aubry  et  al. ,  1988)  or  the  dynamical  evolution 
of  a  jet  (Corke  et  al.,  1994). 

The  aim  of  the  present  study  is  to  apply  such  proce¬ 
dures  on  the  results  of  the  L.E.S.  of  a  3D,  plane  mixing 
layer  spatially  developing  downstream  of  a  flat  plate.  The 
L.E.S.  has  been  performed  on  the  flow  configuration  stud¬ 
ied  experimentally  in  details  in  earlier  works  (see  Delville 
1994,  Ukeiley  1995,  Cordier  1996,  Cordier  et  al.  1997). 
Several  gains  can  be  derived  from  such  a  L.E.S.  computa¬ 
tion:  mainly  the  whole  flow  field  is  known  in  a  deterministic 
way,  thus  the  efficiency  of  the  POD  can  be  analyzed  and 
the  time  evolution  of  the  POD  modes  studied  in  details. 
The  data-base  generated  by  this  simulation  can  be  used 
to  check  assumptions  that  are  required  when  experimental 


procedures  are  applied:  mainly  the  pressure,  experimen¬ 
tally  difficult  to  obtain,  can  be  taken  into  account  and  the 
feedback  relations  (Aubry  et  al.,  1988)  between  the  mod¬ 
eled  structures  and  the  mean  field  can  be  sorted  out. 

We  present  here  the  results  concerning  the  L.E.S.  that 
has  been  performed.  Scalar  “classical”  PODs  and  a 
snapshot-POD  are  derived  from  these  numerical  results. 
It  is  shown  that  the  spatial  eigenfunctions  are  in  good 
agreement  with  the  ones  obtained  experimentally.  Endly  a 
low-order  dynamical  system  is  derived.  The  random  coef¬ 
ficients  obtained  from  the  dynamical  system  and  the  ones 
directly  calculated  from  the  L.E.S.  results  are  compared. 

THE  SNAPSHOT-POD  METHOD 

In  order  to  analyze  the  organization  present  in  the  mix¬ 
ing  layer,  the  Snapshot-POD  method  (Sirovich,  1987)  is  ap¬ 
plied  on  the  numerical  results  given  from  the  L.E.S..  This 
method  based  on  spatial  averaged  quantities,  has  been  pre¬ 
ferred  to  the  “standard”  approach  (Lumley  1967),  since 
simulations  use  an  important  number  of  grid  points  but 
rather  limited  temporal  samples  compared  to  experiments. 

Basic  approach: 

Every  spatio-temporal  event  v{ (x,  t)  is  decomposed  using 
a  mean  and  a  fluctuating  parts: 

=  I7fe) +  »!(£.<)  (1) 

where  the  mean  part  is  expressed  with  an  average  in  time: 

*£(*)=  f  J  V.i{x,t)dt 

The  fluctuating  quantities  can  be  written  using  a  discrete 
basis  of  the  flow: 

npod 

v!i(x,t)=  £  (2) 

n=  1 

where  Npod  is  the  number  of  modes  solved  in  the  POD. 
In  the  snapshot-POD  approach,  the  purely  temporal  eigen- 
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functions  (A^(t))  are  calculated  by  means  of  a  Fredholm 
integral  equation: 

[  C(t,  t ')  il(n)(0  dt’  =  A(n)  A(n\t)  (3) 

Jt 

where  C(t,  tf)  is  the  temporal  velocity-correlations: 

C{t,t')  =  7p  J  2Li(x,t')  dx  (4) 

^  npod 

=  f  £  A(n)(t)^n)(0 

n— 1 

and  where  are  the  (real,  positive)  eigenvalues  of  this 
tensor.  Each  eigenvalue  is  associated  to  the  energy  density 
contained  in  the  corresponding  mode  and  the  sum  of  A^ 
is  equal  to  the  turbulent  kinetic  energy  included  in  the 
integral  domain  ( V ).  The  spatial  eigenfunctions  are  then 
deduced  from: 


The  space-time  coefficients  (a|n)(z,/))  are  then  evaluated 
by  projection: 


£!n)te)  =  j—  f  V t(x,  t)  A<n\t)  dt 


Data  Compression:  “scalar-POD  approach” 

In  order  to  reduce  the  data  storage  requirements,  a  data 
compression  procedure  is  used.  This  compression  is  based 
on  a  scalar  version  of  the  “classical”  POD  introduced  by 
Lumley  (1967).  We  retain  only  “relevant”  modes  of  this  de¬ 
compositions  to  rebuilt  the  flow  realizations  that  are  used 
in  the  snapshot-POD  analysis.  In  this  approach,  we  ex¬ 
press  each  component  of  velocity  of  the  flow  realizations  in 
terms  of  spatial  eigen-functions  depending  only  on  the  x 
and  y  directions: 


-;(->*)  =  f!n)0 x>y) 


*)  =<  Si  >  (*,  V)  +  t)  (7) 

and 

<  V,  >  (*, y)  =  j  J  Vi(x,y,z,t)  dt  dz. 

For  this  scalar  POD,  the  integral  Fredholm  problem  which 
is  solved  is1 

f  f  Raa(x,  y,  x ',  y)  <f^\x',  y’)  dx '  dy  = 

JLX  JLy  “““  (8) 

A(n)  <p£\x,y) 

By  this  way,  instead  of  storing  3  x  Nx  x  Ny  x  Nz  x  Nt 
floating  jxunt  numbers,  we  just  need  to  store  3  x  Nz  x 
Nt  x  N  L  numbers  to  describe  the  temporal  history  of 
the  3D  flow  field.  Here  among  the  N*=NX  x  Ny  modes 
resulting  from  each  of  the  scalar  PODs,  only  NKL  modes 
are  retained.  NKL  is  the  Karhunen-Loeve  dimension  of  the 
scalar  PODs,  i.e.  the  minimum  number  of  modes  necessary 
to  rebuilt  the  turbulent  energy  with  an  accuracy  better 
then  99%. 

In  practice,  the  velocity  correlations  ( Rti(x}  y,  x*,  y‘))  are 
calculated  from  the  L.E.S.  results  using  the  following  rela¬ 
tionship: 

Ra{x,y,x,y')  =  <vi.vi>  (x,y,x',y')  (9) 

-<£,•>  (x,  y).  <  >  (x‘,  y') 

1No  implicit  summation  over  greek  subscripts. 


ain)(M)=/ 

J  Lx  Ly 

and  verify  the  orthogonality  relationship: 


y)  dx  dy  (10) 


Af  [  a<in\z,t).a(im)(z,t) 
Jz *  Jt  l* 


dz  dt  = 


VALIDATION  OF  THE  L.E.S.  RESULTS 
Large  Eddy  Simulation: 

The  L.E.S.  of  the  mixing  layer  is  performed  using  the 
filtered  Navier-Stokes  equations  written  in  the  velocity- 
vorticity  formulation,  by  means  of  the  Reynolds  decompo¬ 
sition  (Lardat  et  al .,  1997).  In  order  to  take  into  account 
the  participation  of  the  small-scale  structures  in  the  fluid 
motion,  the  subgrid  scale  quantities  are  modeled  using  the 
vorticity  transfer  theory  of  Taylor  (1932),  by  means  of  an 
eddy  viscosity  expressed  using  the  mixed-scale  model  (Ta 
Phuoc  1994,  Sagaut  1995,  Lardat  et  al  1997). 

The  spatial  discretization,  based  on  a  finite  difference 
method,  uses  a  M.A.C.  staggered  grid.  The  system  of  equa¬ 
tions  is  solved  in  two  steps  (Lardat  et  al ,  1997).  First  we 
solve  the  transport  equation  of  the  vorticity  by  means  of 
a  Crank- Nicholson,  Adams- Bashforth  approximation  and 
a  QUICK  scheme.  Secondly,  knowing  the  new  vorticity 
field,  the  velocity  is  solved  using  a  projection  (or  fractional 
step)  method,  following  the  work  made  by  Bertagnolio  and 
Daube  (1996).  Let  us  mention  that  this  scheme  is  2nd  order 
accurate  in  both  time  and  space. 

The  L.E.S.  is  performed  on  the  same  flow  configuration 
as  the  experimental  one  obtained  by  Delville  (1995).  The 
velocity  ratio  is  r  =  U2/U1  =  0.59,  where  U1  ~  42.8  m/s 
and  U2  =  25.2  m/s  are  the  magnitudes  of  the  external 
velocities  of  the  boundary  layers  at  the  trailing  edge  of 
the  flat  plate.  The  Reynolds  number,  based  on  and 
on  A U  =  Ui  —  U2i  is  35000  (£Wo  refers  to  the  experimental 
vorticity  thickness  at  a  prescribed  location,  =  30xl0“3 
m).  The  computational  domain  starts  at  the  trailing  edge 
of  the  flat  plate.  It  spreads  over  20  6u0  in  the  streamwise 
direction,  6  £u,0  in  the  inhomogeneous  (vertical)  direction 
(x0  and  lays  over  5  in  the  spanwise  direction  (z). 

The  grid  uses  401x71x55  nodes  in  the  x  x  y  x  z  direc¬ 
tions  and  is  tightened  around  the  centerline  of  the  mixing 
layer  in  the  inhomogeneous  direction  ( y ).  The  flow  is  sup¬ 
posed  to  be  periodic  in  the  spanwise  direction  ( z ).  Using 
this  mesh,  the  grid  filter  width  (A  =  (Az  x  Ay  x  A z)1/s) 
is  close  to  the  Taylor  micro-scale  estimated  following  the 
experimental  results  (Delville,  1995). 

As  regards  the  boundary  conditions  at  the  inflow  surface 
(ie.  the  trailing  edge  of  the  splitting  plate),  the  normal 
component  of  the  mean  velocity  is  initialized  using  two 
turbulent  Whitfield  profiles  for  the  boundary  layers  from 
each  side  of  the  flat  plate.  The  profiles  of  the  tangen¬ 
tial  components  of  the  mean  vorticity  are  then  deduced 
from  the  velocity  profiles.  A  white  noise  is  superimposed 
on  the  tangential  components  of  the  mean  vorticity.  The 
perturbation  magnitude  is  equivalent  to  an  amplitude  of 
10  Ui  on  the  streamwise  velocity  component.  At  the 
outlet  boundary  (rc  =  20  £u;0),  the  tangential  components 
of  the  vorticity  are  calculated  by  extrapolation  along  the 
characteristics.  The  normal  component  of  the  velocity  is 
then  deduced  from  the  vorticity  profiles,  prescribing  that 
the  mass  conservation  is  satisfied.  At  the  upper  and  lower 
surfaces  of  the  domain  ( y  =  ±3  5Wo),  slip  conditions  are 
imposed  (Lardat  et  al,  1997). 
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Figure  1:  Energy  contained  in  the  first  n  modes  of  the 
scalar-PODs. 


To  compare  the  numerical  results  with  the  experimental 
data,  a  portion  of  the  computational  domain  has  been  se¬ 
lected.  This  region  starts  from  5  6o,0  downstream  of  the 
trailing  edge  of  the  flat  plate  and  lays  over  10  £Wo  to  mini¬ 
mize  the  influence  of  the  exit  boundary  condition. 

Validations  of  the  Simulation: 

A  comparison  of  the  numerical  results  has  been  per¬ 
formed  with  the  experimental  data  on  the  mean  and 
Reynolds  stress  profiles.  The  self-similarity  behavior  of 
the  mixing  layer  is  recovered  by  the  L.E.S.  ;  the  vorticity 
thickness  (6^)  and  its  longitudinal  evolution  (dSwjdx)  are 
correctly  predicted  by  the  computation  (L.E.S.:  d6w/dx  = 
4.22  10“2  ;  Experiments:  dS^/dx  =  4.1  10-2).  Very  good 
agreements  are  achieved  on  the  profiles  of  the  Reynolds 
stress  component  (Lardat  et  al .,  1997).  As  regards  the  or¬ 
ganization  of  the  mixing  layer,  the  large  scale  arrangement 
strongly  tridimensionnal,  is  well  predicted.  In  fact,  the  spa¬ 
tial  length  scales  recorded  are  close  to  yx  ~  3.25  8W  in  the 
streamwise  direction  (z)  and  rjz  ~  2/377*  in  the  spanwise 
direction  ( z )  which  are  in  rather  good  agreement  with  the 
ones  generally  admitted  (see  Bernal  and  Roshko,  1986). 


Application  of  the  Scalar-PODs: 

The  scalar-POD  approaches  are  applied  on  the  fluctuat¬ 
ing  velocity  £t  field  provided  by  the  numerical  simulation. 
The  POD  analysis  is  performed  on  1000  temporal  events 
representing  a  dimensionless  time  close  to  T  =  17.  We  ap¬ 
ply  the  POD  on  a  box  with  a  streamwise  extent  of  lOtf^o .  In 
practice,  knowing  an  initial  state,  we  have  first  run  the  cal¬ 
culation  on  1000  time  steps  to  calculate  the  spatial  velocity 
correlations  (Ra(x,  y ,  x\  y'))  (9).  Secondly,  the  eigenvalue- 
eigenfunction  problem  (8)  have  been  solved.  Finally,  the 
simulation  have  been  re-run  from  the  same  previous  ini¬ 
tial  state  to  calculate  the  space-time  coefficients  (0^(2,  <)) 
(10). 

The  convergence  of  the  scalar-PODs  is  presented  on  the 
ratio  A^ /  A^  between  the  energy  contained  in 

the  first  n-POD  modes  and  the  turbulent  kinetic  energy  in 
the  flow  (Fig.  1).  The  convergence  is  relatively  rapid  since 
about  57  %  of  the  turbulent  kinetic  energy  is  contained 
within  the  20  first  modes  and  the  first  520  modes  contain  99 
%  of  the  turbulent  kinetic  energy.  One  can  note  that  a  large 
compression  ratio  (90%)  is  obtained  ;  NxNy/N*  ~  10. 

Iso-contours  of  ^n^(z,y)  are  presented  in  Fig  2  for  four 
selected  modes  (modes  1,  2,  20  and  200)  in  the  plane  (xty). 
In  these  figures,  the  lines  in  bold  correspond  to  positive  val¬ 
ues  while  the  thin  ones  correspond  to  negative  ones.  Con¬ 
sidering  these  eigenfunctions,  one  can  notice  a  streamwise 
shift  when  comparing  modes  1  and  2.  Whatever  the  ve¬ 
locity  components  are,  the  mode  2  is  shifted  downstream 


Figure  2:  Isocontours  of  the  spatial  eigenfunctions 
from  the  scalar-PODs  of  the  three  components  of  the 
velocity  (u,v,w)  for  modes  1,  2,  20  and  200  (from  top 
to  bottom). 


when  compared  to  mode  1.  This  behaviour  is  due  to  the 
convective  nature  of  the  flow  (Rempfer,  1993).  In  fact  the 
first  two  modes  of  the  POD  have  to  be  considered  simul¬ 
taneously  if  one  wants  to  represent  the  mean  convection  in 
the  streamwise  direction.  All  these  eigenfunctions  exhibit 
a  preferred  organization  ;  for  low-order  modes,  this  orga¬ 
nization  follows  the  longitudinal  expansion  of  the  mixing 
layer  which  is  noticeable  on  Fig.  2.  However,  compared 
to  these  low-order  modes,  the  high-order  modes  are  more 
tightened  in  the  vicinity  of  the  centerline  of  the  mixing 
layer  and  correspond  to  smaller  spatial  length  scales. 

The  temporal  evolution  of  the  first  four  random  coeffi¬ 
cients  for  the  scalar  PODs  are  plotted  Fig.  3,  for  the  three 
scalar  functions  ( u ,  v ,  w ),  on  the  middle  z-plane  (2  =  zo). 
A  time  shift  is  visible  between  a^l\zoyt)  and  a^2\zo,t)  and 
between  a(3)(z0,<)  and  a^\zo}t),  which  is  also  due  to  the 
convective  nature  of  the  flow  (Rempfer,  1993). 

The  good  agreement  between  the  computations  and  ex¬ 
periments  can  also  be  shown  when  the  two-point  space  cor¬ 
relations  Ru(x 0,  yoix'^y')  are  considered.  From  the  scalar- 


Figure  3:  Temporal  evolution  of  the  eigenfunctions  of 
the  4  first  modes  of  the  scalar-PODs. 


Figure  4:  Comparison  in  a  plane  (r,  y)  or  (-x,y)  of 
the  experimentally  measured  space  time  correlations 
2/0  (bottom)  with  the  space  correlations  ob¬ 
tained  from  the  contribution  of  the  first  200  modes  of 
the  scalar  PODs  from  the  L.E.S.  (top) 


POD,  this  correlation  can  be  written: 


N* 

Rii(x  o,  yo,  x\  y)  oc  ^2  A(n)£-n)(*o,  0b)£^n)(*',  y) 

n=l 

In  order  to  validate  the  eigenfunctions,  they  can  be  com¬ 
pared  to  the  space-time  correlations  experimentally  deter¬ 
mined  by  using  Taylor  hypothesis  as  in  Delville  (1995).  On 
figure  4,  the  space-time  correlations  experimentally  mea¬ 
sured  are  compared  to  the  contribution  of  the  first  200 
POD  modes  to  the  spatial  correlations.  On  these  figures, 
the  reference  point  (x0)  y0)  is  located  at  the  center  of  the  se¬ 
lected  spatial  domain.  Very  good  agreements  are  achieved. 
We  must  underline  that  the  overall  shape  of  the  measured 
correlations  is  also  well  described  by  considering  only  few 
first  POD  modes.  This  comparison  cross- validates  compu¬ 
tations  as  well  as  experiments.  Note  that  the  u  component 
can  be  described  by  fluctuations  of  opposite  sign  from  part 
to  part  of  the  mixing  layer,  while  the  v  component  can  be 
described  by  fluctuations  in  phase  all  over  the  y  direction, 
alternated  in  sign  in  the  x  direction  (Fig.  4).  Moreover,  the 
shape  of  the  w  component  indicates  the  presence  of  a  strong 
organization  for  the  streamwise  vorticity  component,  this 
organization  being  recovered  from  both  experiments  and 
L.E.S. 

Application  of  the  Snapshot-POD: 

The  snapshot-POD  method  is  then  applied  on  the  fil¬ 
tered  velocity  field  provided  by  the  scalar-PODs  on  the 
first  N  modes.  The  convergence  of  the  snapshot-POD 
is  plotted  on  Fig.  5,  in  the  same  way  as  Fig.  1.  The  conver¬ 
gence  is  also  very  rapid  since  99  %  of  the  turbulent  kinetic 


Figure  5:  Energy  contained  in  the  first  n  modes  of  the 
snapshot  POD. 


Figure  6:  Iso-surface  of  the  spatial  eigenfunctions  of 
the  three  components  of  the  velocity  for  modes  1  (top) 
and  2  (bottom). 


energy  is  contained  within  the  64  first  modes.  However, 
while  the  Karhunen-Loeve  dimension  seems  lower  in  this 
snapshot  POD,  the  relative  number  of  modes  required  to 
contain  99%  of  energy  (NKL/N *=0.064)  remains  of  the 
same  order  of  magnitude  than  in  the  scalar  POD  cases 

(Nkl/N*= 0.1). 

It  is  notable  that  these  eigenvalues  occur  in  pair  of  al¬ 
most  equal  values,  whereas  there  is  a  gap  in  magnitude 
between  these  pairs.  The  eigenvalue  problem  (3)  is  near 
degenerated.  To  understand  the  consequence  of  this  result 
we  present  in  Fig.  6  positive  iso-surfaces  of  $|n)(z,  y,  z) 
for  the  first  two  modes.  Several  general  features  can  be 
noticed  when  considering  these  eigenfunctions.  Whatever 
the  velocity  components  are,  one  can  notice  a  streamwise 
shift  when  comparing  modes  1  and  2.  It  can  be  seen  that, 
essentially,  the  eigenfunctions  of  this  pair  are  representing 
the  same  structure,  one  of  them  just  being  shifted  with 
respect  to  the  other  in  the  streamwise  direction.  The  cor¬ 
responding  expansion  coefficients  A(1)(t)  and  A(2)(<),  are 
found  to  be  analogously  phase  shifted  in  time.  Hence,  if  we 
multiply  each  of  the  eigenfunctions  of  such  a  pair  by  its  cor¬ 
responding  expansion  coefficient,  we  get  a  structure  that 
is  moving  in  the  streamwise  direction.  Analogous  behavior 
is  observed  for  the  first  higher  modes  and  this  behavior  is 
comparable  to  the  one  already  noticed  in  the  scalar-PODs 
spatial  modes  (Fig.  2).  Moreover,  these  eigenfunctions  ex¬ 
hibit  a  preferred  organization  in  the  spanwise  direction. 
Particularly,  $[n\x)  is  clearly  aligned  in  the  spanwise  di¬ 
rection  and  4^,  ^(s)  exhibits  lambda-shape  like  structures 
(Fig.  6).  These  behaviors  are  similar  to  those  observed  in 
the  analysis  of  the  experimental  data  by  Delville  (1994). 
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LOW-ORDER  DYNAMICAL  SYSTEM 

The  goal  sought  after,  is  to  study  the  temporal  behaviour 
of  the  large  scale  coherent  structures.  Therefore,  the  tem¬ 
poral  approach  followed  previously  by  Aubry  et  al  (1988) 
and  Glauser  et  al  (1989)  has  been  employed. 


As  the  pressure  field  is  provided  by  the  numerical  simu¬ 
lations,  a  low-order  dynamical  system  is  developed  by  con¬ 
sidering  the  Navier-Stokes  equations  written  in  pressure- 
velocity  formulation: 

VV  =  0  (11) 

dV  1 

=  -vp+- xrAV 


The  velocity  (U)  and  the  pressure  (P)  are  decomposed 
using  the  Reynolds  decomposition  (1).  The  velocity  fluctu¬ 
ating  part  is  then  written  using  a  discrete  basis  of  the  flow 
following  equation  (2).  This  decomposition  is  introduced  in 
the  transport  equations  (11)  which  are  then  projected  onto 
the  spatial  eigenfunctions  $^(21)  by  means  of  a  Galerkin 
projection.  Since  the  eigenfunctions  are  orthogonal,  a  set 
of  Ngai  time  dependent  ODEs  is  obtained  for  the  so-called 
random  coefficients  of  the  POD  A^(t): 

N9al 

+  P(n)(i) 

m=l 

Ngai  Ngal  V 

+££  Cnmk  [■A(m)(<)j4(k)(i)  -  J4t'T,)(<)J4<*=)(i)j 

171=1  k  —  1 

where  Ngai  is  the  number  of  Galerkin  modes  retained  in 
the  dynamical  system.  The  coefficients  of  the  system  (12) 
bnm ,  Cnmk  are  expressed  as  follows: 

bnm=  (&n\x,y,z),±A&m\x,y,z) 

_VV. $("»)(!,  y,  Z)  -  y,  z).V) 

<W  =  -  (*<»)(.,  y,  z),  x,  y,  z).&k\x,  y,  *)) 

and 

I*n)(t)  =  -  (&n\x,y,z),V?)  =  -  f  ?($(n),n)dS 
'  '  J  QV 

where  n  is  the  external  normal. 

These  ODEs  system  is  solved  using  a  4th  order  Runge- 
Kutta  integration.  At  the  initial  state,  the  random  coeffi¬ 
cients  (A^(t  =  0))  are  initialized  using  the  projection  of 
the  events  j^(x,  t)  at  t  =  0  on  the  spatial  eigenfunctions. 

We  prescribe  A(m)(t)A(k)(t)  =  A^m)  6mk  and  the  bnrn  co- 
efficients  are  estimated  using  the  mean  velocity  given  from 
the  L.E.S..  In  this  first  application,  the  pressure  term  con¬ 
tribution  is  neglected  in  a  fashion  comparable  to  Rajae  et 
al  (1994). 


The  dynamical  system  was  developed  on  the  same  record 
than  the  previous  results.  Two  sets  of  Galerkin  modes 
(Ngai= 10  and  Ngai  =36)  are  retained  to  derive  two  inde¬ 
pendent  dynamical  systems.  For  these  numbers  of  mode, 
the  relative  energy  contained  in  the  system  is  respectively 
42%  and  90%.  The  comparisons  between  the  random  co¬ 
efficients  A^n\t)y  obtained  from  the  resolution  of  the  two 
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Figure  7:  Comparison  of  the  temporal  evolution  of  the 
first  mode  for  the  snapshot  POD  and  the  Dynamical 
Systems  obtained  for  A^a/=10  and  Ngai=36. 


systems  (12),  and  the  projection  of  the  events  u-(x,  t)  onto 
the  spatial  eigenfunctions  $^(x)  are  plotted  for  the  first 
mode  on  the  figure  (7).  When  considering  the  early  evolu¬ 
tion  of  A^x\t)  (t  <  250)  both  dynamical  systems  are  able 
to  follow  the  temporal  evolution  of  this  mode.  The  same 
result  has  been  globally  observed  for  higher  modes.  How¬ 
ever,  if  for  a  short  time  integration,  a  rather  good  agree¬ 
ment  is  achieved,  it  appears  that  for  longer  time,  the  ran¬ 
dom  coefficients  A^n^(t)  obtained  from  the  dynamical  sys¬ 
tems  diverge  from  the  direct  POD  projection.  For  a  small 
truncation  (Ngai= 10),  a  progressive  time  shift  can  be  ob¬ 
served.  When  higher  POD  modes  are  included,  this  time 
shift  can  be  reduced  and  the  typical  period  of  the  modes 
becomes  better  estimated.  Nevertheless,  this  improvement 
is  associated  with  a  divergence  in  the  amplitude  of  the 
first  mode.  This  behaviour  can  be  explained  by  some  of 
the  following  points.  During  the  data  compression  proce¬ 
dure,  obtained  from  the  scalar  PODs,  a  first  truncation  has 
been  performed,  therefore  in  some  sense,  a  certain  amount 
of  small  scales  contribution  has  been  removed  from  the 
L.E.S..  However,  in  projecting  independently  for  each  time 
step  the  “exact”  flow  field  calculated  from  the  L.E.S.  on  the 
spatial  eigenfunctions,  no  bias  is  made  in  the  estimation  of 
the  coefficients  a^n\zyt).  As  regards  the  dynamical  system 
formulation,  neglecting  the  pressure  contribution,  surely 
induces  some  inaccuracy  in  the  model.  In  the  same  way, 
the  instantaneous  contribution  of  the  averaged  field  to  the 
dynamical  system  has  been  simplified  by  considering  only 
its  temporal  average  contribution  =  A^m^m^). 

One  last  reason  for  the  long  term  discrepancy  can  be  re¬ 
lated  to  the  truncation  done  in  the  Galerkin  projection. 
By  taking  into  account  additional  modes  (from  Ngai— 10 
to  36),  still  “energetic”  modes  have  been  added,  and  the 
system  becomes  relatively  less  dissipative.  For  the  dynam¬ 
ical  system  (Ngai= 36),  tests  performed  with  additionnal 
viscosity  have  shown  that  the  time-shift  remains  reduced, 
while  the  magnitudes  of  the  A  coefficients  can  be  lowered. 

CONCLUSIONS 

A  3D  L.E.S.  of  the  plane  mixing  layer  have  been  con¬ 
ducted,  from  which  the  history  of  the  3  velocity  compo¬ 
nents  has  been  calculated  on  a  large  temporal  sample. 
From  the  data  base  generated  by  this  simulation,  two  dif¬ 
ferent  applications  of  POD  have  been  performed.  From  the 
first  POD  approach,  based  on  a  scalar  version  of  the  “clas¬ 
sical”  POD,  the  data  base  has  been  compressed,  allowing  to 
reduce  the  data  storage  requirements  of  90%.  The  spatial 
eigenfunctions  mimic  very  well  the  coherent  structures 
organization  since  a  good  agreement  is  achieved  between 
their  contribution  to  the  space  correlations  and  the  cor¬ 
relations  experiment  ally  obtained.  Using  the  compressed 
data,  a  3  dimensional  snapshot-POD  has  been  performed. 
The  spatial  organization  of  the  eigenfunctions  can  be 
favorably  compared  to  the  experimental  ones. 

Finally,  the  ability  of  the  POD  to  obtain  a  low-order 
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dynamical  system  has  been  underlined.  Two  different  dy¬ 
namical  systems  have  been  performed  using  respectively  10 
POD  modes  and  36  POD  modes.  For  short  time  evolution, 
the  random  coefficients  calculated  from  the  dynamical  sys¬ 
tems  well  agree  with  the  spatial  eigenfunctions  estimated 
from  the  L.E.S.  results.  However,  for  Ngai= 36,  the  dy¬ 
namical  system  are  not  enough  dissipative.  Artificial  dis¬ 
sipation  is  then  necessary  to  follow  correctly  the  long  term 
temporal  evolution  of  the  plane  mixing  layer. 

These  first  results  are  very  encouraging.  In  the  near 
future,  the  influence  of  adding  higher-order  POD  modes 
on  the  flow  structure  representation  will  be  studied  and 
the  role  of  the  pressure  term  clarified. 
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ABSTRACT 

This  numerical  investigation  is  devoted  to  the  direct  sim¬ 
ulation  of  heterogeneous  coaxial  jets.  The  computed  con¬ 
figuration  is  similar  to  the  flow  at  the  exit  of  a  coaxial 
injector  in  the  combustion  chamber  of  the  Vulcain  rocket 
engine  developed  for  Ariane  V.  In  this  flow  pattern,  the 
central  jet  is  composed  of  liquid  oxygen  at  a  low  speed 
and  the  annular  jet  is  constituted  of  gaseous  hydrogen  at  a 
higher  velocity.  This  study  is  led  in  order  to  evaluate  the 
influence  of  the  density  contrast  on  the  pulling  out  of  the 
oxygen  by  hydrogen  for  a  monophase  flow.  The  instabili¬ 
ties  present  in  the  inner  and  outer  shear-layers  involve  the 
presence  of  large  scale  vortices  in  the  near- field.  It  is  found 
that  the  large  scale  two-dimensional  structures  are  at  the 
origin  of  the  destruction  of  the  inner  core.  Its  length  is 
similar  to  this  of  the  liquid  core  in  the  real  two-phase  flow. 
Since  the  breakup  of  the  liquid  jet  is  not  easily  available 
from  experiments  for  this  flow  pattern,  the  simulation  may 
improve  the  knowledge  of  the  atomization  process. 

INTRODUCTION 

Since  several  years  the  simulation  and  the  modelling  of 
turbulent  flows  characterized  by  the  presence  of  organized 
unsteadiness  is  a  research  theme  of  our  team.  In  previous 
studies  (see  Ha  Minh  Sz  Kourta  1993;  Reynier  &;  Ha  Minh 
1996b)  turbulence  models  have  been  applied  to  the  simu¬ 
lation  of  natural  coherent  structures  in  several  flows  such 
as  the  backward-facing  step  and  the  homogeneous  coaxial 
jets.  Here,  the  purpose  is  not  the  simulation  of  coherent 
structures  using  the  semi-deterministic  modelling  (see  Ha 
Minh  Sz  Kourta  1993)  but  their  direct  simulation.  Indeed, 
coaxial  jets  are  characterized  by  the  presence  of  organized 
unsteadiness  in  the  near-field,  already  observed  in  several 
experiments  as  those  of  Gladnick  et  al.  (1990),  Dahm  et 
al.  (1992)  and  Tang  Sz  Ko  (1994).  Coaxial  jets  play  an 
important  role  in  several  industrial  applications  including 
cryotechnic  rocket  engines.  This  flow  pattern  is  character¬ 
ized  by  a  complex  mechanism  of  atomization- vaporization- 
combustion.  As  a  consequence  a  better  understanding  of 
coaxial  twin  fluid  atomizers  is  determinant  to  improve  the 
performance  of  cryotechnic  rocket  engines.  In  this  con¬ 
figuration  the  flow  is  a  two-phase  coaxial  jet  flow  with  a 


central  jet  of  liquid  oxygen  and  an  annular  jet  of  gaseous 
hydrogen.  This  involves  the  presence  of  strong  density  and 
velocity  gradients,  therefore  a  high  shearing  in  the  inner 
mixing  layer  located  at  the  interface  between  the  two  jets. 
The  high  shearing  in  the  near-field  is  at  the  origin  of  the 
presence  of  instability  waves.  Generally  the  destruction  of 
the  liquid  core  is  due  to  the  generation  of  instabilities  that 
initially  distort  and  finally  break  the  core  as  observed  ex¬ 
perimentally  by  Andrews  (1993),  Issac  et  al.  (1994)  and 
Mayer  (1994).  For  Andrews  (1993),  the  perturbations  may 
be  generated  by  two  mechanisms:  an  helical  wave  owing  to 
the  relaxation  of  the  velocity  profile  within  the  core  and  the 
coherent  self-induction  of  larger  perturbations  from  smaller 
ones  by  vortex  chaining  and  nonlinear  interaction.  Accord¬ 
ing  to  Lefebvre  (1992)  in  liquid  of  low  viscosity,  instability 
waves  protrude  into  the  coflowing  stream.  These  waves 
enhance  the  interaction  between  the  gas  and  the  liquid, 
causing  the  crests  of  the  waves  to  detach  and  form  liga¬ 
ments,  which  break  down  into  drops.  As  a  consequence 
the  size  of  the  drops  is  directly  related  to  the  diameter 
of  the  ligaments.  Other  studies  as  those  of  Tsai  Sz  Viers 
(1992),  Wu  et  al.  (1992)  and  Wu  &  Faeth  (1995)  are  in 
agreement  with  this  process  of  drop  formation.  These  au¬ 
thors  reported  that  the  formation  of  the  drops  is  associated 
to  the  turbulent  eddies.  This  process  is  important  because 
it  determines  the  initial  properties  of  the  dispersed  phase, 
which  influences  mixing  rates,  secondary  breakup,  which 
represents  the  conversion  of  large  drops  into  droplets,  col¬ 
lisions  and  separated  flow  within  the  dispersed  region  (e.g. 
Wu  et  al.  1992).  The  mechanism  of  droplet  formation  is 
very  difficult  to  study  in  cryotechnic  rocket  engines.  Ac¬ 
cording  to  Andrews  (1993)  this  is  due  to  the  optical  density 
associated  to  the  atomization  process.  The  investigation  of 
Ingebo  (1992)  on  the  breakup  of  a  cryogenic  liquid  jet  re¬ 
ports  that  the  mechanism  of  breakup  is  much  more  difficult 
to  study  than  the  atomization  of  fuel  or  water  jets,  because 
the  surface  temperatures  were  often  near  the  boiling  point 
of  the  liquid. 

In  the  perspective  of  a  better  understanding  of  the 
breakup  of  the  liquid  core,  the  present  paper  focuses  on 
the  simulation  of  two-dimensional  structures  in  heteroge¬ 
neous  coaxial  jets  for  a  monophase  flow.  This  approach 
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of  a  two-phase  flow  can  amaze  at  first  sight,  however  it 
is  valid  under  certain  conditions.  This  way  is  acceptable 
if  the  effects  of  the  surface  tension  forces  are  negligible 
against  the  aerodynamic  effects.  In  the  injection  of  pro¬ 
pellants  at  a  supercritical  pressure,  surface  tension  forces 
play  a  less  important  role  than  the  aerodynamic  effects. 
Moreover,  this  assumption  is  valid  at  lower  pressures  when 
high  gradients  of  density  and  velocity  are  present.  This 
is  confirmed  by  the  experiments  of  Gicquel  et  al.  (1995) 
and  Lefebvre  (1992)  which  showed  that  the  atomization  is 
highly  dependent  on  both  momentum  and  gas-liquid  den¬ 
sity  ratio.  In  fact  the  flows  where  this  approach  is  ac¬ 
ceptable  are  the  flows  at  high  Weber  numbers.  This  way 
has  been  applied  already  to  the  simulation  of  heteroge¬ 
neous  coaxial  jets  using  turbulence  models  by  Reynier  & 
Ha  Minh  (1996a).  These  previous  simulations  have  allowed 
a  good  prediction  of  the  length  of  the  liquid  core  which  is 
according  to  Wu  et  al.  (1992)  of  the  same  order  that  the 
length  of  the  potential  core  in  an  monophase  flow. 

GOVERNING  EQUATIONS 

This  numerical  study  focuses  on  the  pulling  out  of  an 
heavy  fluid  by  a  lighter  one  under  the  action  of  the  large 
scale  structures  evolving  in  the  flow.  In  coaxial  jets,  three- 
dimensional  effects  are  weak  up  to  the  end  of  the  potential 
core.  According  to  Tang  &  Ko  (1994)  in  air- air  coaxial 
jets,  azimuthal  structures  involved  by  the  bifurcation  of 
two-dimensional  vortices  of  the  near  region  are  not  present 
before  the  intermediate  region  located  downstream  the  end 
of  the  inner  potential  core.  As  the  objective  of  this  work 
is  the  approach  of  the  breakup  of  a  liquid  jet,  the  study 
focuses  on  the  near-field  of  the  flow.  Therefore,  a  two- 
dimensional  approach  is  used  for  the  calculations.  The 
problem  deals  with  coaxial  jets,  so  the  Navier-Stokes  equa¬ 
tions  and  the  equation  for  the  concentration  rate  of  oxygen 
are  solved  in  axisymmetric  coordinates.  These  equations 
are  given  as: 

Equation  for  the  conservation  of  the  density  p: 
d  d  i  Q 

-p+—pu  +  ~rpv  =  0  (1) 

Equation  for  the  conservation  of  the  momentum  pU : 


Figure  1;  Computed  configuration  of  hydrogen-oxygen 
coaxial  jets 


Equation  for  the  conservation  of  the  total  energy  E: 

State  equation  for  a  perfect  gas 

P  =  (y-  1  )pei  (9) 

Equation  for  the  conservation  of  the  mass  fraction  rate  of 
oxygen  pA: 


§ifA  +  i!;{l’AV--k§;A} 

Where  Pr  and  Sc  are  respectively 
Schmidt  numbers. 


(10) 

the  Prandtl  and 


dtpU+9x  [pUl  +  p 

+7 -o 

Equation  for  the  conservation  of  the  momentum  pV : 


hpV  +  h{puv  ~  pSxT] 

+  rlrr{pV’2  +  P_/t5rr}  =  7{i’-'‘5e9} 

(3) 

The  terms  of  the  strain  tensor  are  given  by  the  following 
expressions: 

5  -odJL  2  (dU  ,dV  V\ 
brr~2dr  ~3  [d^  +  -d7  +  7) 

(4) 

c  0dU  2  fdU  ,  dV  V\ 

5~  =  2dT-3br  +  i>r  +  7) 

(5) 

5  -5  _9U  dV 

Sxr-Srx-—  +  — 

(6) 

See=2V_2(8U  +  aV  +  V) 
r  3  \  dx  dr  r  ) 

(7) 

NUMERICAL  TOOLS 

The  finite- volume  method  proposed  by  MacCormack 
(1981)  is  used  for  the  simulations.  It  is  an  implicit-explicit 
scheme,  but  to  simulate  the  natural  unsteadiness  of  coaxial 
jets  only  the  explicit  part  is  used  in  this  study.  The  nu¬ 
merical  scheme  uses  the  predictor-corrector  step  procedure 
and  the  equations  are  solved  under  a  conservative  form. 
The  method  is  accurate  to  the  second  order  in  time  and 
space. 

The  flow  to  be  computed  is  presented  in  Figure  1.  It  is 
an  heterogeneous  flow  of  oxygen  and  hydrogen  discharging 
in  a  hydrogen  atmosphere.  The  density  ratio  between  the 
two  flows  is  equal  to  16.25.  The  diameter  D  of  the  injector 
is  equal  to  8.76  mm.  The  diameter  of  the  central  jet  is 
5.14  mm.  The  velocity  Ui  in  the  round  jet  is  10  m.s”\  the 
Reynolds^number  is  3617  and  the  density  p0 2  is  equal  to 
1.3  kg.m  3.  The  inner  diameter  of  the  annular  jet  is  6.74 
mm  and  the  outer  diameter  is  8.76  mm.  The  lip  thick¬ 
ness  between  the  two  jets  is  0.5  mm.  The  velocity  U Q  in 
the  annular  jet  is  160  m.s  1  and  the  Reynolds  number  is 
3880  for  a  density  of  hydrogen  pH 2  equal  to  0.08  kg.m”3. 
The  exit  temperature  is  300°K,  the  pressure  Pc  at  the  noz¬ 
zle  and  in  the  computational  field  is  initially  equal  to  one 
atmosphere. 

The  computational  domain  extends  over  22  diameters  in 
the  streamwise  direction  and  6.8  diameters  in  the  radial 
direction.  The  grid  is  158  x  97.  Along  x,  the  mesh  is  thin 
and  uniform  with  140  cells  up  to  12  diameters  then  coarser 
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up  to  22D.  In  the  radial  direction  the  grid  uses  a  stretching 
outside  the  flow. 

At  the  nozzle,  the  velocity  profile  has  been  calculated 
with  the  same  numerical  code  which  has  been  adapted  to 
the  simulation  of  a  flow  in  a  pipe.  In  the  rest  of  the  inlet 
section  a  wall  is  present  (see  Figure  1)  then  homogeneous 
Dirichlet  conditions  are  applied  to  velocity  and  homoge¬ 
neous  Neumann  conditions  are  imposed  for  concentration 
rate,  pressure  and  density.  As  the  upper  boundary  is  far 
from  the  flow  homogeneous  Neumann  conditions  are  ap¬ 
plied  for  all  quantities.  The  lower  boundary  is  located  on 
the  flow  axis,  then  symmetry  conditions  are  applied.  The 
outlet  conditions  are  deduced  from  the  characteristic  re¬ 
lationships.  These  boundary  conditions  are  non-reflective 
and  have  been  originally  developed  by  Thompson  (1987) 
for  the  resolution  of  the  Euler  equations. 

RESULTS 

In  Figure  2  and  3  the  profiles  of  the  streamwise  and  radial 
velocities  are  presented  for  four  sections  of  the  flow  located 
at  0.25D,  ID,  3D  and  8D  for  the  streamwise  velocity  and 
0.25D,  ID,  2D  and  8D  for  the  radial  velocity.  The  results 
have  been  averaged  over  110000  time  steps.  The  profiles 
of  the  streamwise  velocity  show  a  strong  decrease  of  the 
velocity  in  the  annular  flow  near  the  nozzle.  At  0.25D  the 
velocity  is  equal  to  45  per  cent  of  the  exit  velocity.  This  de¬ 
crease  can  be  explained  by  the  strong  contrasts  of  density 
and  velocity  between  the  two  flows.  As  the  hydrogen  flow 
carries  filaments  of  oxygen,  there  is  a  high  transfer  of  mo¬ 
mentum  from  the  hydrogen  to  the  oxygen.  This  is  support 
by  the  profiles  of  the  radial  velocity  (see  Figure  3)  which 
establish  the  existence  of  a  high  level  of  radial  velocity  at 
0.25D.  Downstream,  the  profiles  of  the  longitudinal  veloc¬ 
ity  become  flat  and  at  8D  the  maximum  is  located  on  the 
axis.  The  time-averaged  profiles  of  the  spanwise  velocity 
put  in  evidence  a  strong  decrease  of  this  quantity  between 
0.25D  and  ID.  Next  this  location  the  time- aver  aged  level 
of  the  radial  velocity  is  weak. 

The  unsteady  field  of  the  mass  fraction  rate  of  oxygen 
is  represented  in  Figures  4  to  9.  This  figures  show  only 
the  near-field  of  the  flow  and  the  corresponding  domain 
extends  over  12  diameters  along  x  and  3  diameters  along 
y.  The  times  corresponding  to  this  figures  are  11.75  ms 
for  the  figure  4,  12.4  ms  for  the  figure  5,  13.06  ms  for  the 
figure  6,  13.71  ms  for  the  figure  7,  14.36  ms  for  the  figure 
8  and  15.02  ms  for  the  figure  9.  The  visualizations  put  in 
evidence  a  high  unsteadiness  of  this  flow  which  is  involved 
by  the  shearing  occuring  in  the  inner  mixing  layer  located 
between  the  two  jets.  The  shearing  involves  the  presence  of 
large  scale  vortices  attested  by  the  high  value  of  the  time- 
averaged  radial  velocity  (see  Figure  3)  in  the  near-field.  A 
previous  numerical  study  (see  Reynier  &  Ha  Minh  1996a) 
has  established  the  coupling  between  the  large  scale  vor¬ 
tices  occuring  in  the  inner  mixing  layer  and  the  pulling  out 
of  the  oxygen.  This  coupling  between  the  pulling  out  of 
oxygen  and  the  presence  of  coherent  structures  is  in  agree¬ 
ment  with  the  atomization  model  proposed  by  Andrews 
(1993)  where  the  destruction  of  the  liquid  core  is  involved 
by  the  presence  of  coherent  structures  at  the  liquid-gas 
interface.  Here,  the  flow  has  not  a  high  degree  of  orga¬ 
nization  notably  after  3D.  This  lack  of  organization  can 
be  explained  by  the  high  contrast  of  density  between  the 
two  streams.  The  high  level  of  instability  involves  a  fast 
destruction  of  the  core  of  oxygen.  As  the  flow  moves  down¬ 
stream  the  mixing  becomes  more  and  more  intense  particu¬ 
larly  after  the  destruction  of  the  oxygen  core  at  approxima¬ 
tive^  2D.  This  high  mixing  involves  the  pulling  out  of  the 
oxygen  in  the  shape  of  filaments  then  the  tearing  of  these 
filaments.  The  mechanism  of  destruction  of  the  oxygen 
core  is  involved  by  the  vortices  dominating  the  inner  mix¬ 
ing  layer.  They  involve  the  formation  of  filaments  but  also 
the  tearing  of  oxygen  packets  along  x  on  the  axis.  The  de¬ 


Figure  2:  Time-averaged  profiles  of  streamwise  veloc¬ 


ity 


y/D 

Figure  3:  Time-averaged  profiles  of  radial  velocity 


tachment  of  the  oxygen  packets  occurs  at  approximatively 
2D,  then  they  are  destroyed  as  they  move  downstream. 
This  phenomenon  of  destruction  of  the  oxygen  core  is  very 
violent.  This  agrees  qualitatively  with  the  experiment  of 
Rupe  (1962)  who  observed  that  laminar  flows  can  be  more 
unstable  and  are  atomized  in  a  more  violent  fashion  than 
fully-developed  turbulent  jets,  due  to  higher  radial  velocity. 

CONCLUSION 

In  a  previous  numerical  study  (see  Reynier  &  Ha  Minh 
1996a)  turbulence  models  have  been  applied  to  the  sim¬ 
ulation  of  natural  coherent  unsteadiness  in  heterogeneous 
coaxial  jets.  These  previous  results  showed  that  the  large 
scale  structures  evolving  in  the  flow  are  at  the  source  of  the 
pulling  out  of  the  oxygen  by  the  hydrogen  stream.  Here, 
the  direct  simulation  of  heterogeneous  jets  confirms  this 
point.  The  simulation  of  the  destruction  process  of  the 
oxygen  jet  is  very  violent.  This  result  is  in  qualitative 
agreement  with  past  experiments.  To  precise  this  process 
and  to  approach  the  real  flow  which  is  a  two-phase  one, 
with  a  better  reliability,  three-dimensional  calculations  will 
be  necessary. 
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ABSTRACT 

This  paper  presents  an  experimental  study  on  the  interaction  of 
a  heated  jet  with  a  normal  shock  wave  in  a  Mach  1 .6  nozzle.  The 
jet  is  heated  to  maintain  a  uniform  Mach  number  ahead  of  the 
shock  wave.  The  jump  in  the  TKE  across  the  shock  is  studied  by 
means  of  2D-LDV  measurements.  Turbulence  in  the  central  part 
of  the  jet  is  amplified.  Amplification  factors  are  close  to  those 
predicted  by  the  linear  theory  and  do  not  depend  on  the  jet 
temperature.  In  the  mixing  layer,  the  jump  in  the  turbulence 
properties  are  strongly  dependent  on  the  jet  temperature.  This 
seems  to  be  primarily  due  to  changes  in  the  mean  shear  across 
the  shock.  These  results  are  discussed  in  the  light  of  available 
theoretical  results  in  the  field. 

INTRODUCTION 

The  problem  of  turbulence-shock  wave  interaction  is  an 
important  topic  for  transonic  and  supersonic  aeronautics.  Two 
most  important  basic  types  of  flows  in  the  domain  are 
shock/boundary  layer  interactions  and  under-expanded 
supersonic  jets.  Only  shock/boundary  layer  interactions  have 
been  widely  studied  (see  Delery  &  Marvin,  1986,  for  a  review). 

In  the  past  ten  years,  a  renewal  of  interest  in  this  field  has  been 
motivated  by  projects  relating  to  high  capacity  supersonic 
transportation.  The  problem  has  been  revisited  on  a  more 
fundamental  basis.  In  particular,  the  case  of  homogeneous 
turbulence,  which  is  propitious  to  theoretical,  numerical  and 
experimental  developments,  has  been  investigated.  Experiments, 
direct  numerical  simulations  and  linear  theory  have  brought  a 
deeper  comprehension  of  the  way  a  shock  wave  transforms 
incoming  homogeneous  turbulence  (see  e.g.  Lele,  1994  for  a 
review).  A  unified  phenomenology  is  not  yet  available,  but 
several  elements  of  the  puzzle  are  now  understood  : 

i)  Results  of  turbulence  amplification  by  a  shock  wave  are 
mediated  by  compressible  mechanisms  which  correspond  to 
rapid  exchanges  between  potential  and  kinetic  energy  in  the 
nearfield  behind  the  shock  surface.  Experimental  investigation 
of  these  mechanisms  by  conventional  means  (LDV  and  hot  wire) 
is  a  difficult  task.  They  are  described  by  the  "linear  interaction 
analysis"  (LI A)  developed  by  Ribner  in  the  early  fifties  (see 
Ribner,  1953)  and  have  been  well  characterised  by  direct 
numerical  simulation  (see  Lele,  1994,  for  a  review). 

ii)  Compressible  mechanisms  make  the  interaction  between 
solenoidal  isotropic  homogeneous  turbulence  and  a  shock  wave 
less  energetic  than  a  rapid  one-dimensional  compression  of 
similar  strength  (see  Lee  et  al.,  1991,  Jacquin  et  al.,  1993(a)). 

iii)  As  shown  by  linear  theory  and  DNS,  the  issue  of  a 
turbulence/shock  interaction  depends  on  anisotropy  of  the 
upstream  turbulence  (see  Mahesh  et  al.,  1993)  and  on  the 
thermodynamic  properties  of  the  upstream  fluctuating  fields 


(Ribner,  1953,  Moore,  1954,  Hannapel  &  Friedrich,  1995, 
Mahesh  et  al.  1995,  1996).  Concerning  the  last  problem,  recent 
theoretical  and  numerical  results  have  shown  that  upstream 
entropy  fluctuations  may  have  a  large  influence  on  the 
interaction  (Mahesh  et  al.,  1997).  This  will  be  especially 
addressed  in  the  present  paper. 

The  question  is  whether  these  results  will  help  to  explain 
complex  flows.  For  instance,  in  shock/boundary  layer 
interaction,  amplification  of  turbulent  kinetic  energy  (TKE)  is 
usually  stronger  than  that  predicted  by  linear  analysis  of 
homogeneous  turbulence.  Several  mechanisms  can  promote 
turbulence  in  such  a  flow.  First,  the  highest  turbulence 
amplifications  reported  concern  cases  where  the  boundary  layer 
separates  in  the  shock  region.  In  this  case,  turbulence 
amplification  is  mainly  due  to  separation.  This  is  not  a  direct 
effect  of  the  shock  itself.  Moreover,  when  penetrating  into  the 
boundary  layer,  the  shock  induced  compression  is  spread  within 
a  lambda  shaped  region.  The  resulting  bulk  compression  may 
lead  to  higher  amplification  than  that  obtain  behind  a  pure  shock 
wave  (see  ii)).  The  anisotropy  of  the  incoming  turbulence  and 
shock  unsteadiness  are  other  possible  sources  of  amplification. 
Also,  as  shown  by  Mahesh  et  al.  (1997),  the  thermodynamic 
equilibrium  which  prevails  in  an  adiabatic  boundary  layer  may 
promote  amplification  through  baroclinic  effects. 

A  new  contribution  to  the  field  is  presented  here.  It  concerns 
the  case  of  a  heated  jet  interacting  with  a  normal  shock  wave. 
This  experiment  aims  at  providing  new  results  on  the  topic  of 
free  turbulent  flow/shock  interaction.  Several  arguments  evoked 
above  in  the  case  of  shock/boundary  layer  interaction  will  be 
discussed  in  the  framework  of  this  new  study. 

EXPERIMENTAL  SET-UP 

The  experimental  set-up  is  shown  in  figure  1.  A  cylindrical 
sonic  jet  discharges  into  the  sonic  section  of  a  M0=  1.6 
rectangular  nozzle  whose  test  section  is  120x100mm2  (subscript 
0  denotes  the  freestream,  and  j  the  jet  flow).  The  jet  exit 
diameter  is  Dj  The  exit  static  pressure  of  the  jet  is 

matched  with  that  of  the  freestream.  Consequently,  the  jet  and 
the  freestream  expand  together  from  Mj  =  Mext  =  1  to  Mach 

1.6.  Due  to  the  change  in  the  nozzle  aspect  ratio  between  the 
sonic  throat  and  the  test  section,  the  jet  becomes  slightly 
elliptical.  The  Mach  number  being  uniform  across  the  jet,  a 
pressure  adapted  interaction  between  the  jet  and  a  shock  wave 
can  be  produced. 

The  velocity  difference  between  the  jet  and  the  external  flow  is 
modulated  by  changing  the  jet  temperature.  An  electric  heater, 
of  diameter  33mm,  located  upstream  of  the  ejector,  is  used  for 
this  purpose.  Reduction  in  the  diameter  from  this  value  to 
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Dj  -  17 mm  is  made  by  means  of  a  converging  ejector  (see 

figure  1).  Three  values  of  the  total  temperature  ratio  have  been 
tested  RT=Tij/Tie„=  1,2,3  (2^=300*). 

A  normal  shock,  controlled  by  a  second  throat,  is  placed  in  the 
diverging  part  of  the  nozzle  at  location  xc  =  110mm.  This 
corresponds  to  xcJDj  =6.5  .  Figure  2  shows  a  schlieren  picture 

of  the  flow  for  RT  = 3  (7*,-  =900 K)1. 

MEAN  FIELD 

Figure  3  shows  the  evolution  of  the  mean  velocity  along  the  jet 
axis  for  RT  =3  .  The  flow  accelerates  in  the  diverging  part  of  the 

nozzle  (from  M}=\  at  x=0  to  Afj=1.6  at  jr~90mm  ).  Then  the 
velocity  of  the  jet  decreases  progressively  as  it  develops.  When 
the  interaction  with  the  shock  wave  occurs,  the  jet  still  possesses 
a  «  potential  »  core.  Figure  3  also  shows  that  just  after  the  shock 
deceleration,  the  flow  undergoes  a  slight  reacceleration.  This 
comes  from  the  separation  of  the  nozzle  boundary  layers.  This 
phenomenon  is  visible  on  the  Schlieren  picture  in  figure  2. 

Figure  4  shows  the  Mach  number  and  the  total  temperature 
distributions.  Measurements  were  obtained  by  means  of  pressure 
and  thermocouple  probes  situated  10mm  upstream  (x=  100mm) 
and  downstream  (x=  120mm)  of  the  shock.  The  upstream  part  of 
the  flow  is  characterized  by  a  uniform  Mach  number 
distribution,  except  for  a  slight  deficit  in  the  mixing  layer  region 
due  to  the  pressure  loss  in  the  jet  generator  boundary  layer.  The 
uniformity  of  the  pressure  and  Mach  number  distributions  are 
preserved  across  the  shock.  The  photograph  in  figure  2  shows 
that  the  shock  surface  is  not  distorted  by  the  jet. 

In  the  present  flow,  there  are  variations  of  the  mean  shear  in 
the  diverging  part  of  the  nozzle  as  well  as  through  the  shock. 
This  will  influence  the  property  of  turbulence  inside  the  mixing 
layer.  These  variations  may  be  evaluated  by  considering  the 
vorticity  equation  : 

dc o  -  — _  _  1  - 

*—  =  0).  vu-u>divu+^-'Vpx  Vp  (1) 

dt  p2 

in  the  case  of  a  one  unidirectional  compressible  shear  flow,  such 
that : 

jcoKO.nc^.o),  u  =  (u(x,z),o,o) 

|p  =  p{x\  p  =  p (x,z),  M  =  M(x) 


By  combining  (1)  and  (2)  and  expressing  the  baroclinic  term  by 
means  of  Euler  and  state  equations,  the  solution  of  (1)  can  be 
easily  written  as  : 


d__ 

dx 


=  0 


(3) 


Consequently,  the  mean  shear  in  the  jet  mixing  layer,  £2  = - , 

dz 

increases  in  the  diverging  part  of  the  nozzle  and  decreases 
through  the  shock  wave2.  If  one  discards  entropy  effects,  these 
variations  only  depend  on  the  Mach  number  distribution.  From 
the  sonic  section  to  the  shock  location  (where  Af-1.6),  it  is 
found  that  the  mean  shear  increases  by  a  factor  1.8  (this  value  is 
only  indicative  because  the  spreading  of  the  mixing  layer  is  not 
accounted  for  in  this  model).  Note  that  the  cylindrical  shear 


Turbulence  visible  outside  the  jet  in  the  downstream  part  of  the  flow 
corresponds  to  the  boundary  layers  that  develop  on  the  side  windows. 

Note  that  in  a  one-dimensional  compression  of  homogeneous  shear 
flow,  one  has  £IU  =  Cte  (see  Mahesh  et  al.  1994).  This  corresponds  to 
an  inverse  evolution  with  respect  to  equation  (3).  Consequently, 
homogeneous  linear  theory  is  not  very  useful  for  the  interpretation  of 
the  jet  mixing  layer  behavior  in  the  present  experiment. 


layer  which  hits  the  shock  wave  is  not  in  equilibrium.  It 
corresponds  to  the  initial  development  region  of  a  jet  and, 
moreover,  it  is  accelerated  by  the  nozzle.  Through  the  shock 
wave,  one  obtains : 


£l=£l=C(Mo)= 

n0 


2  +  (y-l)  Ml 
(Y  +  l  )M% 


(4) 


where  subscript  0  denotes  flow  conditions  ahead  of  the  shock, 
and  subscript  1,  those  behind.  ForM0=1.6,  Qj/Qq  =  0.49  :  the 

mean  shear  is  reduced  by  nearly  a  factor  2.  One  can  expect  that 
this  change  in  the  mean  shear  across  the  shock  will  lead  to  a 
reduction  of  turbulence  in  the  downstream  part  of  the  flow. 

Finally,  since  the  Mach  number  is  uniform  across  the  flow,  the 
velocity  ratio  Ry^U^/Uj  ,  the  density  ratio  =Pj/pext  and 


the  convective  Mach  number  Mc  of  the  jet  mixing  layer  only 


depend  on  the  total  temperature  ratio  Rr  =  Ti{  IT  and  on  the 

1  l  1  ext 

Mach  number : 


Rjj  —  Rt  1  2 ,  Rp  =  Rj  1 ,  Mc — M — — 


1/2 


-1 


V/2+i 


(1) 


with  M  -  Mj  =  Mext .  When  M=  1 .6,  this  gives  : 


Rj 

7?p 

Mc 

l 

l 

1 

0 

2 

0,5 

0,71 

0,27 

3 

0.33 

0.58 

0.42 

Table  1 :  Mixing  parameters  for  M-L6 

Mc  does  not  exceed  0.5  so  that  compressibility  induced  effects 
in  the  mixing  layer  are  assumed  to  be  weak. 

SECOND  ORDER  STATISTICS 

Mean  and  fluctuating  velocity  components  have  been  measured 
using  a  two-component  LDV  system.  MgO  particles  were  used 

for  seeding  the  flow.  The  Reynolds  stresses  u'2  ,  w'2  and  mV 
(w*  is  the  vertical  component  of  velocity)  have  been  determined 
for  the  three  jet  temperatures,  with  and  without  shock  wave. 
Vertical  profiles  measured  ahead  (*=100mm)  and  downstream 
(*=110 mm)  of  the  shock  wave  are  drawn  in  figure  5.  Figure  6 

shows  the  corresponding  amplification  factors  R  ?  =  m’?/m^, 

u  w  u 

RW2  ~  w'i  and  /^2=<Zi2Ao  where  q2-  u'2+lw'2  . 

In  figure  5,  it  is  observed  that  even  though  the  mixing  layer  has 
not  reached  the  jet  axis  at  *=100  mm  ( x/Dj  =  6.5 ),  the  central 

region  of  the  jet  is  turbulent.  This  turbulence  is  created  by  the 
heating  system  located  inside  the  jet  generator  and  is  locally 
quasi-homogeneous.  The  flow  amounts  to  the  superposition  of  a 
turbulent  mixing  layer  and  a  quasi-homogeneous  turbulence 
located  in  its  centre.  Figures  5  and  6  show  that  the  issue  of  the 
shock  interaction  depends  upon  the  region  considered. 

Quasi  homogeneous  turbulence 

In  the  jet  axis  region  ( -5 mm  <  z  <  5 mm  )  turbulence  is 
amplified.  Values  of  the  energy  components  ,  the  anisotropy 

ratio  ( ),  the  upstream  turbulence  rate  (<?o/^o  ) 

corresponding  amplification  factors  obtained  at  z  =  0  are 
reported  in  table  2. 
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u'1  is  amplified  whereas  wu  is  weakly  affected.  It  is 
remarkable  to  note  the  invariance  of  R  2  =  wf  /m*^  with  respect 

to  Rt  .  The  value  obtained  is  very  close  to  that  given  by  the 
linear  theory  (LIA)  (for  initial  solenoidal  isotropic  turbulence). 


Rj 

V2 

-o'2 

4 

Vo 

t  1 

1  is. 

wf 

4 

4 

IL 

(mV2) 

(mV2' 

)4 

)  Vo 

“0 

ofo 

4 

& 

1 

166 

235 

0.71 

455 

0.055 

1.67 

1.0 

1.18 

1.17 

2 

433 

450 

0.96 

640 

0.057 

1.72 

1.0 

1.65 

1.23 

3 

625 

590 

1.06 

770 

0.055 

1.67 

1.1 

1.50 

1.35 

Table  2  :  Turbulence  properties  on  the  jet  axis  ahead  (subscript  0, 
x  =100  mm)  and  behind  (subscript  1,  x-120  mm)  the  shock. 


However,  LIA  predicts  a  slight  increase  in  w'2  which  is  not 
observed.  Variation  in  the  TKE  amplification  factor 

R  2  =  Q\  I  %  mainly  reflects  changes  in  the  initial  anisotropy 
ratio  m'q/ Wq  3.  It  may  be  noted  in  table  2  that  the  initial 

qa 

turbulence  rate  — 3 ~  is  constant.  The  upstream  TKE  value  q0 
U  o 

thus  increases  by  a  factor  -VT  from  7?^  =1  to  3.  This  does  not 

seem  to  have  a  significant  influence. 

On  the  whole,  one  may  conclude  that  in  the  central  region  of 
the  jet,  quasi-homogeneous  turbulence  is  subjected  to  a 
moderate  amplification  and  that  the  level  of  this  amplification  is 
compatible  with  linear  theory.  It  does  not  seem  to  depend 
neither  on  the  total  temperature  of  the  flow  nor  on  the  initial 
turbulence  level. 

Mixing  layer 

As  shown  in  figures  5  and  6,  changes  in  the  turbulent 
quantities  in  the  mixing  layer  are  multiform. 

For  Rt  =  1 ,  turbulence  in  the  wake  region  almost  behaves  like 

2  2 

that  in  the  central  region  :  u'  is  amplified,  whereas  W  is 
almost  unaffected,  as  well  as  u'w ' .  In  figure  6,  there  is  no 
significant  change  in  the  amplification  factors  for  Rr  =  1 , 

except  in  the  external  region  of  the  jet  ( \z\  >  1 1  mm )  where  they 

increase.  An  increase  in  the  external  region  is  also  observed  for 
the  two  other  temperature  cases.  It  is  provoked  by  a  slight 
thickening  of  the  jet  behind  the  shock  combined  with  the  low 
levels  of  initial  turbulence  that  prevails  in  this  region. 

For  RT  *  1 ,  changes  in  the  mixing  layer  are  characterized  by  a 

decrease  of  the  shear  stress  m’w’  ,  an  increase  of  w'2  and  by  a 


3  As  explain  above,  the  turbulence  in  the  jet's  central  region  is  produced 
inside  the  jet  generator.  It  is  subjected  to  an  axisymmetric  compression 
inside  the  sonic  ejector  and  to  an  axisymmetric  expansion  inside  the 
tunnel  nozzle.  The  first  effect  leads  to  a  strong  decrease  of  the 

anisotropic  ratio  u'  /  W  ,  the  second,  to  an  increase.  The  values  of 

u'qJ w'o  obtained  at  the  shock  level  results  from  the  combination  of 
these  mechanisms.  This  depends  on  the  jet  velocity,  that  is  on  the  jet 
temperature. 


spectacular  reduction  of  u'2  .  Interpretation  of  these  results  is 
not  straightforward.  Several  mechanisms  must  be  accounted  for. 
The  first  one  is  the  reduction  in  the  mean  shear  across  the  shock. 
The  second  is  possible  temperature  effects.  The  third  concerns 
possible  measurement  bias  due  to  LDV  particle  dynamics. 

As  shown  above,  the  mean  shear  is  reduced  by  a  factor  of  two 
across  the  shock.  This  means  that  turbulence  production  is 
reduced  in  the  same  proportion.  The  shock  wave  transforms  the 
supersonic  mixing  layer  into  a  subsonic  one  which  is  less 
energetic  but  which  contains  intense  turbulent  fluctuations 
which  results  from  the  shock  interaction.  Consequently,  a  region 
exists  downstream  of  the  shock  where  this  turbulence  must  be 
dissipated  in  order  to  match  new  equilibrium  conditions. 

One  can  try  to  characterize  the  above  mechanisms  by 
considering  the  normalised  Reynolds  stresses 

Ou/AU  ,  Ow/A U  ,  Ou/ow  and  iCw'/q2  where  a  denotes  a 

r  m.s  (e.g.  au  =  ft  )  and  AU  y  the  difference  in  the  mean 
velocity  across  the  mixing  layer  (AU - Uax js  - £/ext ;  AU  varies 
according  to  (3) :  AU/U  =  Cte).  Figure  7  shows  the  evolution 

of  these  quantities  along  a  line  connecting  points  where  q2  is 
maximum.  We  focus  on  the  region  70mm  <  x(mm)  <  150mm 
between  the  nozzle  acceleration  and  the  development  region  of 
the  jet. 

In  the  supersonic  regime,  ahead  of  the  shock,  one  obtains 
Ou/AU  ~  0.2,  aw/AU  =0.1  and  Tw'/q1  =0.19  for  RT=  3.  A 
slightly  higher  value  of  aw/AU  is  obtained  for  RT  =2.  For  the 
anisotropy  ratio,  one  obtains  ou/gw  =1.9  for  RT  =3  and  1.7  for 
Rt=  2.  Experiments  on  plane  supersonic  mixing  layers  at 
moderate  convective  Mach  numbers  suggest  smaller  values  for 
GU/AU  (==0.175),  and  similar  values  for  aw/AU  (see  a 

review  in  Smits  &  Dussauge,  1996).  Consequently,  the 
anisotropy  ratio  of  the  supersonic  upstream  in  the  present 
experiment  are  higher  than  that  commonly  observed. 

The  normalised  components  of  energy  increases  behind  the 
shock4  and  the  anisotropy  is  decreased  to  a  more  standard  value 
of  1.3.  It  is  striking  to  note  that,  compared  to  results  on 
equilibrium  plane  mixing  layer,  the  flow  at  x=150mm  is 
characterised  by  higher  values  of  the  normalised  components  of 
energy  but  by  a  more  conventional  value  of  the  anisotropy  ratio. 
As  discussed  by  Jacquin  et  al.  1995,  high  values  of  the 
anisotropy  ratio  have  been  already  observed  in  LDV 
measurements  in  supersonic  heated  jets,  but  also  in  plane  mixing 
layers.  As  explained  by  these  authors,  seeding  bias  could 
contribute  to  this  behaviour.  If  any,  these  bias  should  depend  on 
the  mean  shear  amplitude.  This  could  explain  why  the 
anisotropy  ratio  increases  upstream  with  RT  and  why  it 

decrease  suddenly  in  the  subsonic  regime  behind  the  shock.  This 
problem  is  the  object  of  ongoing  research. 

As  seen  through  the  behaviour  of  the  normalised  TKE 

production  ( u'w'/q 2  )  adaptation  of  the  flow  to  downstream 

conditions  takes  place  at  a  distance  of  20  or  30mm.  A 
characteristic  distance  for  the  dissipation  can  be  evaluated  as 


cu/AU  is  sujected  to  a  sharp  overshoot  at  the  shock  level.  This 

overshoot  is  commonly  observed  in  LDV  results.  It  is  due  to  a 
combination  of  the  shock  intermittency  and  particle  drag  effects 
(differences  in  the  dynamic  response  of  particles  of  various  sizes  to  the 
shock  induced  deceleration).  The  amplitude  and  the  thickness  of  this 
spurious  "particle  drag  turbulence"  is  enhanced  by  shock  unsteadiness 
(see  Jacquin  et  al.  1993(b)  for  further  details). 
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(Ax)e  =  U  te  where  U  denotes  the  convection  velocity  and 

Te-(#2/2)/e  the  dissipation  time  scale,  with  e  the  TKE 

dissipation  rate.  The  time  scale  ( Te )  can  be  only  evaluated 
ahead  of  the  shock  if  equilibrium  between  production  and 
dissipation  is  assumed.  The  dissipation  rate  (Bq  )  then 

reads:  e0  =  (u'w'dU/dz^.  The  change  in  the  time  scale 

=  (*2/2)/e  through  the  shock  is  unknown.  Assuming  that 

Teis  conserved5,  that  is  Tel  =  ze0  =  using  the 

maximum  values  of  mV  and  dU/dz  at  x=100 mm  and  the 
velocity  U\  measured  downstream  in  the  mixing  layer,  one 
obtains  (Ax)£l  =  UxTz0^60mm  for  RT  =  2  and  (Ax)£l  ~  40mm 
for  Rt  =  3 ,  respectively.  These  values  are  equivalent  to  those 
observed.  This  confirms  that  turbulence  of  the  mixing  layer  is 
subjected  to  rapid  dissipation  behind  the  shock.  However,  as 
shown  in  figure  7,  properties  of  the  downstream  mixing  layer 
strongly  differ  from  those  upstream.  This  behavior  is  not  yet 
understood. 

ABOUT  TEMPERATURE  EFFECTS 
As  previously  stated,  a  theoretical  analysis  of  the  behaviour  of 
the  fluctuating  field  across  the  shock  may  be  done  in  the 
framework  of  Ribner  s  linear  interaction  analysis  (LIA).  Ribner 
(1953)  developed  the  case  of  upstream  pure  shear  waves.  This 
analysis  has  been  recently  extended  by  Mahesh  et  al.  (1996)  to 
cases  where  additional  entropy  waves  are  present  (see  also 
Mahesh  et  al.,  1997).  In  the  latter  case,  the  theory  shows  that  the 
correlation  between  the  vorticity  and  the  entropy  fluctuations 
strongly  influences  the  evolution  of  turbulence  across  the  shock. 
In  particular,  Mahesh  et  al.  have  shown  that  a  negative 
temperature- velocity  correlation  m'7”  upstream  of  the  shock 

enhances  TKE  amplification;  whereas  if  mT  is  positive  the 
amplification  is  strongly  damped. 

Considering  the  homogeneous  region  of  the  flow,  one  may 
assume  that  correlation  between  velocity  and  temperature  is  not 
promoted  by  increase  in  the  mean  total  temperature  of  the  flow. 
This  is  consistent  with  an  invariance  of  the  amplification  factors 
with  respect  to  RT  as  observed  in  figure  6. 

Contrarily,  temperature  effects  may  be  present  in  the  mixing 
layer.  The  mean  temperature  gradient  across  the  mixing  layer  is 
aligned  with  the  mean  velocity  gradient,  so  that  strong  positive 
u'T  must  prevail6.  Thus  baroclinic  effects  could  contribute  to 
lessen  the  TKE  amplification  in  the  mixing  layer.  Unfortunately, 
this  effect  cannot  be  isolated  from  the  consequences  of  the  mean 
shear  reduction.  One  can  however  note  in  figure  7  a  slight 
reduction  in  the  downstream  TKE  as  Rr  increases. 


Actually,  DNS  and  experiments  on  homogeneous  turbulence  show  that 
the  shock  compression  leads  an  augmention  of  enstrophy  (see  work  by 
Lele  et  col.  and  Jacquin  et  al.,  1993(b)).  This  must  lead  to  an  increase  of 

6  A  fluid  particle  displaced  impulsively  from  the  free  flow  into  the  jet 
flow  possesses  a  negative  velocity  fluctuation,  u\  and  a  negative 
temperature  fluctuation  T;  both  fluctuations  are  positive  in  the  opposite 
case. 


SUMMARY 

An  experiment  is  presented  in  which  a  heated  jet  interacted 
with  a  normal  shock  wave.  A  pressure  matched  interaction  with 
a  nearly  uniform  Mach  number  distribution  ahead  and  behind 
the  shock  wave  was  obtained  by  heating  the  jet.  The  turbulent 
properties  of  this  interaction  were  measured  by  LDV.  The 
following  trends  have  been  observed  : 

1)  The  present  flow  is  characterised  by  two  types  of  turbulence. 
The  first  is  a  quasi-homogeneous  turbulence  convected  in  the 
central  region  of  the  jet.  The  second  is  that  which  develops  in 
the  mixing  layer. 

2)  The  quasi-homogeneous  turbulence  prevailing  in  the  central 
part  of  the  flow  is  amplified  at  a  rate  which  corresponds  to  the 
linear  theory  prediction.  Amplification  is  not  significantly 
affected  by  temperature  nor  by  the  upstream  turbulence 
intensity. 

3)  The  reduction  in  the  mean  shear  which  is  induced  by  the 
shock  wave  leads  to  a  decrease  in  the  turbulence  activity  of  the 
subsonic  mixing  layer  downstream  of  the  shock.  Accordingly,  a 
reduction  in  the  TKE  may  be  found  through  the  shock  in  the 
mixing  layer. 

4)  Individual  components  of  the  TKE  measured  upstream  of 
the  shock  are  characterized  by  a  significant  departure  from 
equilibrium.  Equilibrium  is  never  reached  downstream  of  the 
shock,  perhaps  because  of  the  jet  development. 

5)  As  suggested  in  Mahesh  et  al.  (1997),  baroclinic  effects 
could  participate  in  the  reduction  of  the  TKE  amplification  in 
the  shear  layer.  Unfortunately,  this  cannot  be  corroborated  in  the 
present  experiment  because  the  flow  behind  the  shock  is 
dominated  by  rapid  changes  in  the  TKE  production  and  by  a 
rapid  dissipation. 

6)  On  the  whole,  this  jet/shock  interaction  is  characterized  by 
TKE  amplifications  smaller  than  those  observed  in  other  flows 
where  shock-induced  distortions  are  present,  as  in 
shock/boundary  layer  interactions  or  in  under-expanded  jets. 
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Figure  1 :  Experimental  set-up 


Figure  2  :  Schlieren  picture  (20  ns)  of  the  flow 


- * -  x- 100  mm 

- • -  la  130mm 


- e -  x  ■  100  mm 

- • -  xa  120mm 


Figure  4  :  Mach  number  and  total  temperature 

distribution  at  x=  100mm  and  120  mm  for  the  case 

Rt  =  3 


(a) 


Figure  3  :  Evolution  of  the  mean  velocity  along  the  jet  axis  for  the  case 
%  =  3. 
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